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the polymers pass within Angstroms of each other leads to a decrease of the electronic energy 
of the combined system that is significantly larger than the thermal energy and thus promotes 
interchain aggregation. We consider the competition of this attractive interaction with the screened 
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I. INTRODUCTION 

Conjugated polymers have been the subject of intense 
research in the physics, chemistry, and materials science 
communities. These molecules provide both a laboratory 
to explore one-dimensional conductors Q, |^ |^ and the 
potential for a plethora of applications in such disparate 
devices as solar cells j^] and chemo/biosensors^ J^X 8J. 
The simplest conducting polymer system is polyacety- 
lene. The underlying cause of its conductivity, however, 
is somewhat subtle owing to the fact that, due to the 
Peierls instability!^ 0| in one dimension, this molecule 
appears to be a semiconductor (based on an analysis of 
uncorrelated electrons). Here we consider the case of an 
electron-doped polymer in which the Fermi energy lies 
within the conduction band. 

The utility of conjugated polymers is generally real- 
ized in such complex chemical systems created by adding 
functional groups to the basic, conjugated backbone. 
Such modified polymers typically include side groups to 
modify the conductivity of the backbone by doping in 
either electrons or holes; polar or ionizable side groups 
can also dramatically increase the solubility of these 
molecules in water Fluorescent, water-solubilized 

conjugated polymers such as polyparaphenylene vinylene 
(PPV) hold promise for use as biosensors since their 
fluorescence is readily quenched by electron acceptors 
such as methyl viologen (MV). Biomolecules tagged with 
MV can be detected optically at remarkably low concen- 
tration using these water soluble conjugated polymers. 
Fluorescence quenching can be quantified by the Stern- 
Volmer constant, -f^svi which describes the loss of quan- 
tum efficiency per mole of quencher present. Since the 
Ksv - lO^M-i for the PPV/MV system, MV tagged 
biomolecules can be detected via fluorescence quenching 
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at 100 nanomolar concentrations'l2]. Before this analytic 
tool can be realized, a number of complications must be 
addressed. These include understanding and accounting 
for the dramatic effect of "bystander" molecules such as 
surfactants on the value of ^s\/E3|, as well as address- 
ing the role played by chain length of the PPV [13 . [Till . 
Additionally, the PPV/MV system is hindered as an an- 
alytical tool by the tendency of the PPV to aggregate 
in solution and to bind nonspecifically with other small 
molecules. Neutron and light scattering^^ demonstrate 
that PPV solutions form large aggregates under low-salt 
conditions. With added salt, the scattering data can be 
fit by rod-like structures having a persistence length of 
80nm 10 suggesting the bundling of the chains in solu- 
tion. 

Understanding the tendency of metallic polymers to 
aggregate in solution is the central focus of this paper. 
We propose an interchain adhesion mechanism specific to 
conjugated polymers that should apply to doped, solubi- 
lizcd PPV and other conducting polymers. When two 
polymer chains approach each other to Angstrom-scale 
distances at some point along their backbones, the inter- 
chain tunnelling of the electrons (holes) in the conduction 
bands of the two molecules decreases the total electronic 
energy of the system. This decrease in the energy of the 
electronic degrees of freedom is primarily due to the cre- 
ation of a low energy, localized state at the crossing point 
of the two chains. These states lead to the aggregation of 
the two chains since, for physically reasonable parame- 
ters the energy decrease per crossing point is on the order 
of a few /cbT. Such a electron tunnelling mechanism has 
been previously discussed as a source of an intrachain, at- 
tractive interaction leading to polymer collapse in good 
solvent 

The mechanism that we propose can be thought of as 
analogous to the creation of a molecular covalent bond 
with the distinction that in the present system, the bond- 
ing and anti-bonding states that are localized at the 
crossing point of the two chains are not created from 
normal atomic energy eigenstates, but rather are pulled 
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out of the delocalized states of the charge carriers on the 
chains. In the remainder of the article, we first explore 
the consequences of this type of interchain bonding (in 
section in three stages of increasing complexity. We 
consider the effect of a single chain-crossing point that 
leads to interchain tunnelling in section III Al We then 
turn to the case of multiple crossing points by consider- 
ing an ordered array of such points in section lll Bl before 
turning to the more physical case of a thermally disor- 
dered set of crossing points in 111 CI Finally, we consider 
the attraction between two polymers aligned in a parallel 
configuration in IIIDI 

In order to determine the adhesive properties of these 
conjugated polymers, we need to compare the decrease 
in the energy of the electronic degrees of freedom of the 
molecule, to the increase of chain conformational free en- 
ergy associated with adhesion. We turn to the latter cal- 
culation in section imi before combining these two parts 
of the analysis in section IIVI wherein we discuss the re- 
sults and propose experimental tests of this work. 

II. INTERCHAIN ADHESION 

Since we suggest that the adhesive properties of the 
conducting polymers are a generic consequence of pop- 
ulating the conduction band of the these extended 
molecules, we develop our theory based on a minimalistic, 
tight-binding Hamiltonian for these conduction electrons 
on a chain of N sites of the form 

N/2 

Ho = -t ^ K>(^+1| + 1^+1X^1, (1) 
e=-N/2 

where |^) is a state vector for an electron on the tight 
binding site. In the absence of strong electron-electron 
correlations!^ such a single particle approach is justified 
and many properties of such conjugated polymers have 
been predicted via such simple tight- binding models |l9j . 
The overlap integral t is not precisely known for a num- 
ber of these systems, but recent spectroscopic results on 
chemically related systems 19] suggest that the t is on 
the order of 1 — AeV. 

To discuss the modification of the electronic states of 
this tight-binding model due to the proximity of another 
such polymer at one point or a set of points {ii}, we 
modify Eq. ^ by introducing another copy and including 
an interaction term between the chains parameterized by 
a second hopping matrix element t': 

N 

^0 = 5] 5](K+1,J)(^,J| + K,J)(^ + 1,J|)(2) 

i=l,2fcl 

Hj = -t' ^ {\ll){l2\ + \l2){ll\). (3) 

In the above equations, the sum on j is over the two 
chains, while the interaction Hamiltonian, Hj allows for 



interchain hopping at a selected subset of sites along the 
polymer, {£/}. The interchain hopping matrix element 
has been shown numerically to be on the order of O.leV 
for similar chemical systems '2^. It is essential to note 
that since the interchain tunnelling matrix element is ex- 
ponentially sensitive to the interchain separation, it is ac- 
ceptable to suppose that this matrix element will be finite 
at only isolated points. In fact, we will later show in sec- 
tion lllll that the positions of these "crossing-points" to be 
spatially correlated or anticorrelated along the arc length 
of the chain due to chain conformational free energy. A 
simple pictorial representation of the basic problem for 
the case of only one hopping site is shown in figure^ 




FIG. 1: A pictorial representation of the simplest crossed- 
chains configuration in which the inter-chain hopping is al- 
lowed at the central site, £ = 0. The spheres represent 
the tight-binding sites for the conduction electrons along the 
chains. The chains need not be straight as shown here. 

Because of the large separation of time scales be- 
tween the electronic and conformational degrees of free- 
dom of the polymer chain, we may ignore conforma- 
tional changes in the polymer while discussing the mod- 
ification of the electronic states due to the proximity of 
the two chains at a given point or set of points. To 
show that we may ignore the thermal motion of the 
point of closest approach in discussing the electronic 
structure of the molecules, we compare the time scale 
for monomer motion over distances of a Bohr radius 
(the distance over which one expects t' to vary rapidly) 
Tchain o^B \/ TTi / {ksT) to characteristic times for elec- 
tronic structure reorganization Tdocton fi/t. We find 
that Toioctron/Tchain ~ 10"^'^ SO the adiabatic approxi- 
mation is eminently reasonable. 

Similarly, we note that the temperature of the sys- 
tem (typically 3QQK or less) is significantly less than the 
Fermi temperature of electronic system so in the remain- 
der of the calculation we will determine the total free 
energy of the electronic degrees of the freedom in a zero 
temperature approximation whereas we will consider the 
role of temperature and consequently entropy while dis- 
cussing the conformational and translational degrees of 
freedom of the polymer. We now develop our calculation 
for the binding energy of the two chains by computing 
separately the decrease in the energy of the electronic de- 
grees of freedom due to interchain tunnelling and increase 
in chain free energy due to the loss of some translational 
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and conformational entropy. From the sum of these two 
effects, we determine the effective binding energy. First 
we consider a single crossing point as shown in figure Q 

A. A Single Crossing Point 

We take the single crossing point to lie in the middle 
of both chains {ij} — {0} as shown in figure H We re- 
turn to crossing configurations of lower symmetry later. 
We note that the interaction Hamiltonian, Eq. |31 can be 
diagonalized in the basis of symmetrized and antisym- 
metrized chain occupation states: 

\i,±) = ^[\i,l)±\i,2)]. (4) 

so that in this basis Eq. |21 takes the form 

Hi = -t' E (5) 

It is immediately clear that the action of the Hamilto- 
nian on the subspaces of the antisymmetrized and sym- 
metrized states is identical except for the exchange of 
t' — > —t'. We discuss the energies of states in both 
subspaces in parallel. The Hamiltonian is also symmet- 
ric under the parity operator, P\£,zL) = \ — £,±); since 
only the states symmetric under P will be affected by 
the crossing point, we focus on these even parity states 
in the following. We make an ansatz for the unnormal- 
ized, (anti-) symmetrized, even parity eigenstates of the 
Hamiltonian by writing 

m 

|fc,±>- E co8{ak\£\+cbf)\£) (6) 

£— — rn 

where (p^ represents the wavevector-dependent break in 
the phase of the wavefunction due to the interchain in- 
teraction on chains of = 2m + 1 tight binding sites. 
Boundary conditions require that the amplitude of the 
wavefunction at sites ±{111+1) vanishes, which with Eq.^l 
leads to 

afc(m + l) + 0^^ ^(Y^\ peZ. (7) 

One can show that the states \k, ±) with arbitrary phase 
(j)^ satisfy the time-independent Schrodinger equation all 
sites away from the crossing point since, via direct calcu- 
lation, one finds 

E^{£^0,±\k,±) = {£^0,±\H\k,±), (8) 

with 

= -2tcos{ak). (9) 

However, to simultaneously satisfy this eigenvalue equa- 
tion at the crossing point as well as Eq.Q we must choose 
the phase angle 0^ such that 

t' 

sin(fca) = =F— tan [fca(m + 1)] . (10) 



The solutions of this transcendental equation that deter- 
mine the allowed wavevectors k (and energies via Eq. 
are shown graphically as the set of intersections in fig- 
ure 121 where the left and right hand sides of Eq. ^| are 
plotted. 




ka 

FIG. 2: Graphical solution of Eq. 1101 Allowed k values 
lie at the intersections of the solid line with the long dashes 
(symmetric band) and with the short dashes (anti-symmetric 
band). The vertical lines represent the unperturbed k values. 

The principal consequence of the crossing point is the 
creation of two localized bound states out of the sym- 
metrized/anitsymmetrized states ±). We examine the 
symmetrized states first. The conditions on the tun- 
nelling matrix element t' for the appearance of this bound 
state can be inferred from Eq. ^| If 

1<M(™ + 1) (11) 

then there is no solution for real wavevector ka < 
7r/[2(m + 1)], but a new solution appears along the imag- 
inary axis in the complex k plane at fc = in, where k is 
given by 

sinh(aK) = — tanh [aK{m + 1)] . (12) 

From this relationship and Eq. |51 we determine the en- 
ergy of the resulting bound state to be 



Eb = -2*^' 1 + ^ tanh2[K(m + 1)] (13) 

where k is the solution from Eq.^l From the magnitude 
of the imaginary wavevector, it is clear that this bound 
state is localized on the length scale of at/t' ^ a. It 
is important to note that, due to the appearance of this 
bound state on any chain of reasonable length, the energy 
of the electronic degrees of freedom is decreased by a 
quantity on the order of the tunnelling matrix element 
t' . In the limit of an infinite chain where m — > 00, one 
sees from the above equations that an ~ 0{1) and the 
tanh(-) in Ea. ll3l becomes one in agreement with previous 
work |21|. This rearrangement of the electronic degrees 
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of freedom has dramatic consequences for the polymer 
conformational dynamics as we show below. 

In addition to this lower energy bound state, a second 
localized state is pulled from the conduction band of an- 
tisymmetrized states. In this case t' — > —t' as discussed 
below Eq.El From Eas. llllO we now find that the com- 
plex wavevector associated with the bound state becomes 
k = ZK + tt; the bound state appears at the edge of the 
Brillouin zone and the energy of the state is —Ef,. It is 
interesting to note that the appearance of these two lo- 
calized states at energies ±Et, is directly analogous to the 
appearance of bonding/anitbonding states in covalently 
bonded atoms^- In the present case, however, we build 
these states out of spatially extended conduction-band 
states rather than the localized atomic orbitals. 

There are still other consequences of interchain tun- 
nelling for the scattering states that remain in the con- 
duction band. In other words, all the even symmetry, 
extended states of the unperturbed chains are shifted in 
energy due to their interaction with the crossing point 
at site zero. These scattering states remain extended, 
i.e. retain a purely real wavevector, but that wavevec- 
tor shifts due to the interaction with the crossing point. 
For t' = we trivially find these extended states at 



k 



(0) _ 



7r(2p -I- 1)/{2N) for integer p. Due to scatter- 



ing at the crossing point, these wavevectors shift so that 
kp — > kp + Akp. These shifts can be computed by ex- 
panding both sides of Eq.^|near fep"-* . From this expan- 
sion we find to leading order in t' /N that 



Afc. 



2tN sin kp 



1 



t' cos kp 
2tN sin^ kr 



-O 
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Using the above shifts in kp, we can compute using Eq.|U| 
the sum of the shifts in the energy levels of these scat- 
tering states in the symmetrized band. Recalling that 
we need to also include the analogous shifts in the chain 
antisymmetrized band, which eliminates terms odd in t', 
we find that the total shift is 



AE=-[- 
N 



2 (mt) , 

cos kp 



2t sin^ fc„ ' 



(15) 



where the sum is over all filled energy levels from the 
bottom of the conduction band to rrif set by the Fermi 
level of the system. Since the sum of all such energy shifts 
of the extended states vanishes in the limit of large N as 
1 /N, the changes in the energy of the scattering states in 
the band is not significant for long chains. Hereafter we 
will ignore the finite length {0{1/N)) The principal result 
is that there is an attractive interchain interaction due 
to a reorganization of the electronic degrees of freedom 
of the system. The energetically dominant part of this 
reorganization is the appearance of the "bonding" bound 
state that lowers the electronic energy essentially by t''^/t. 

Since the dominant contribution to the electronic en- 
ergy shifts comes from a localized state at the crossing 
point, it is perhaps not surprising, that, when the above 



analysis is generalized to the case of a crossing point at an 
arbitrary point along the chain, the principal, qualitative 
result is unchanged. We will show below that this ba- 
sic result is effectively insensitive to interactions between 
such crossing-points in chains that form a multicrossing- 
point "braids". We do this in two steps by first con- 
sidering an ordered array of crossing points and then by 
examining a spatially disordered set of crossing points-see 
figure 121 



B. An Ordered Array of Crossing Points 

One may wonder about the effect of a multiple set of 
interchain crossings. Specifically, if, as it appears above, 
one crossing point decreases the overall energy of the elec- 
tronic states of the system, is it more profitable for the 
chains to form many such crossing points, and if so, at 
what density? To address this question, we first con- 
sider the most simple such structure - an ordered array 
of crossing points at every M'^ site along the chains. The 
ordered structure forms, in effect, a superlattice in which 
the unit cell is constructed from one crossing point and 
a basis of M — 1 tight binding sites between that con- 
secutive crossing points. Qualitatively, one expects the 
localized states centered at each crossing point to merge 
into what may be termed an impurity band in the lattice. 
We explore this point below paying particular attention 
to the interchain binding energy per unit length of the 
polymers. 




•C)aX)CXXX)CX3»(XXX»C)CKX»»(X)CXDCXX)CXX« 




•ooocxxccoocaD^ooooa^^ooooococoooc)* 

'S. ^ 




FIG. 3: Two polymers interacting at a series of localized sites. 
In the upper figure the two chains form an ordered structure 
in which a crossing point appears after exactly M tight bind- 
ing sites. The lower figure represents the disordered version 
of the braid. Below each braid is a pictorial representation of 
the appropriate tight-binding Hamiltonian in which the filled 
circles represent the interchain tunnelling sites. In the disor- 
dered two-chain braid, the number of sites between crossing 
points varies along the chains, M — > Mi, but the same Mi 
sites exist on both chains between the crossing points. 
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To begin we note that within a single supercell of the 
superlattice, the tight-binding Hamiltonian becomes 



AI~1 



Hn - ^ {K,n)(^+l,n| + |£ + l,n)(^,n|} 
e=i 

+X\M,n){M,n\ (16) 

where state \£, n) is a position eigenstate representing an 
electron at the site in the n}^ supercell. Due to the 
interchain tunnelling at the M"^ site in each supercell, 
there is a localized potential for the electron to reside 
there. In the above equation and in what follows we 
find it convenient to suppress the ± notation describing 
the interchain symmetrized and antisymmetrized states. 
We have seen above in Eq. |S1 that the only diS'erence in 
the form of the Hamiltonian acting on these two sub- 
spaces of opposite chain-exchange symmetry is the shift: 
t' — > —t'. In the above equation, that t' dependent term 
in the Hamiltonian is given in terms of A which takes on 
the value t' (— i') for the antisymmetrized (symmetrized) 
states. To complete our superlattice Hamiltonian, we in- 
clude a simple tunnelling term to couple the states in each 
supercell. Here, as before, we assume that the presence of 
the crossing chain does not affect the intrachain hopping 
matrix element, although such effects could be taken into 
account using the present tight-binding formalism. The 
coupling term is 



-t{|Af,n)(l,n+ 1| + \l,n+l)(M,n\}. (17) 



The first term in the above equation couples the crossing- 
point site in the n"^ supercell to the tight-binding site to 
its right (in the (n -|- 1)"^ supercell) with the usual in- 
trachain hopping matrix element, while the second term 
allows for electron hopping from the first site in {n+ 1)*'' 
supercell onto the crossing point site to its immediate left 
that is part of the n}^ supercell. The full Hamiltonian of 
the system is simply the sum of these two terms summed 
over all N supercells: 



N 



(18) 



where we have used periodic boundary conditions so that 
the {N + 1)* supercell is identified with the first super- 
cell. Since principal effect of the interchain tunnelling on 
the energetics of the system has been shown to be the 
appearance of localized bound states for a single crossing 
point, the net effect of periodic boundary conditions used 
here is minimal. 

In order to diagonalize the Hamiltonian given by 
Eq. El we introduce two wavevectors conjugate to the 
real-space degrees of freedom that index the supercell 
and the site within each supercell by writing unnormal- 
ized states 



N M 



\q, k) = Y.Y. exp^'"*'" (e"" + A.fee"'''"') 1^- 



n=l e=i 



where Aq,k will be determined below. The allowed values 
of q are fixed by the periodic boundary conditions to be 
q = 2Trn/NMa where n is an integer such that: —N/2 < 
n < N/2. We now insist that these states, which are 
eigenstates of the supercell translation operator are also 
energy eigenstates: H\q, k) — E{q, k)\q, k). This relation 
can be translated into three related equations formed by 
projecting the above relation onto three carefully chosen 
states. We require that 

{ln\H\q,k) = E{q,k){ln\q,k) (20) 
{l,n\H\q,k) = E{q,k){ln\q,k) (21) 
{M,n\H\q,k) - E{q,k)(J.,n\q,k), (22) 

where i ^ 1, Af and the index of the supercell n is arbi- 
trary. The first condition determines the energy values 
to be E{q,k) — —2tcos{ka). Note that since condition 
Eq. 1201 applies to the interior of each supercell, the en- 
ergy eigenvalues appear to be independent of q and A. 
This apparent lack of q- and A-dependence is false, as 
we will soon be forced to require that the allowed val- 
ues of k depend on these parameters. The second and 
third conditions above acknowledge the special role of 
the crossing-point in the Hamiltonian. From the second 
condition, we determine the parameter Aq^k- lu order 
to satisfy the eigenvector/eigenvalue equation at the left- 
most site in a supercell, we find 



^i(kaM-q) _ ^ 



I _ g-i{kaM+q) 



(23) 



In the above result we have used the energy eigenvalues 
obtained from enforcing the first condition. Finally, the 
combination of the third condition, Eq.|221along with the 
energy eigenvalue used above leads to a relation between, 
k, q, and A of the form 



— tan [kMa) 



1 



cos{qMa) 
cos{kMa) 



sin(fca). (24) 



The above equation may be read as an implicit form of 
the relation of k upon q and A, i.e. k = k{q,X). Note, 
that in the limit of a very sparse array of crossing-points, 
i.e. M — > oo we recover the condition on imaginary k 
for the isolated bound state given in Eq.^l However, as 
M becomes smaller, we find a small modification in this 
condition and the introduction of a lattice momentum 
(q) dependence in k. We now turn to the calculation 
of the density of states g{E) = dAf/dE where J\f is the 
total number of energy eigenstates with energy less than 
E. Using the density we compute the total energy of the 
electronic degrees of freedom of the braided polymers. 
The density of states can be written as 



1 



dk If^ 



(25) 



(19) 



The second derivative in the product above is computed 
trivially from the dispersion relation. Using the quan- 
tization of the wavevector q, the first derivative can be 
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expressed as: 

dJ\f dN dn dq 



dk 



dn dq dk 



= M ■ 



NMa dq 
2tt ' dk 



point is highly independent of the local density of such 
crossing points at least until M decreases to order unity. 



(26) 



where the first term in the product reflects that fact that 
for each wavevector (? = (?„ = 2TTn/{Na) there are M al- 
lowed values of fc, one for each of the M bands in the sys- 
tem. The second term in the above product expresses the 
g-quantization condition; this leaves only the last term to 
be computed from Eq. 1241 This derivative is given by 



1.5 



g(E) 



0.5 



dq 
dk 



2t 



-Mcos{kaM) + sin(fcaM)£2£(^(fca) 



M am{kaM) 



CPA 
Ordered 



My 1 - (57 sin(fcaM) csc(fca) -I- cos(fcaM)) 

(27) 

Now, the total energy of the electronic states of the chains 
can be computed from Eqs.|5Sl|5Hl and[57| 

One expects from the calculation of the isolated 
crossing-point that the change in electronic energy is al- 
most entirely due to the creation of a narrow band of 
bonding states below the conduction band so it is rea- 
sonable to imagine that the total change in electronic en- 
ergy is independent of the Fermi level. For convenience 
we choose the Fermi level to be Ep = 0. Since there is at 
most one energy eigenstate for each (complex) value of fc, 
we may write the total energy of the electronic system as 
an integral over k of the density of states and the disper- 
sion relation. All possible states below the Fermi level 
are explored by integrating along the following contour 
in the complex plane: C — (— zoo, 0) lj(0, 7r/2) which is 
the combination of the negative imaginary axis and the 
segment of the real axis spanning to 7r/2. The density 
of states along this contour is shown in figure01as a func- 
tion of energy. Lastly, using Eq.[^wc note that for there 
to be an energy eigenstate with a given {q, k), it must be 
true that |A/(2i) sin(MaA:) csc(afc) -I- cos(A:aAf)| < 1 so 
that the denominator of Eg. 1271 must be real. To enforce 
this condition over the k integration we take the real part 
of the integral to write 



E= Re [ ^E(k)dk. 
Jc ok 



(28) 



The change in the total electronic energy as a func- 
tion of t' and M can then be computed; these results are 
shown in figures [S] and The observed even-odd effect 
in figure m is an artifact of our choice of the Fermi level. 
For chains with even M, a band gap will open at Ep = 
resulting in an additional reduction in the energy. In 
practice, we expect that the Fermi level of doped poly- 
mers will not lie at a place of such high symmetry and the 
dominant contribution to the energy will be the creation 
of localized bound states. 

The principal result of this calculation is that, due to 
the highly localized nature of the bonding states at each 
crossing-point, the total interchain binding energy is es- 
sentially a linear function of the crossing point density. 
In other words the attractive interaction at each crossing 



■ 0.5 
E/2t 



0.5 



FIG. 4: The density of states as a function of energy for an 
ordered array of crossing points (solid line) and a disordered 
array of crossing points (dashed line). In both cases the den- 
sity of binding sites is p = 1/M = 0.125 and \/{2i) = -0.25. 



C. A Disordered Array of Crossing Points 

As a final point regarding the electronic contribution 
to the interchain adhesion, we now allow the set of cross- 
ing points to become disordered in a specific manner as 
shown in the lower part of figure|21 We allow the number 
of tight-binding sites between two crossing points (M) 
now to vary along the chain in an uncorrelated manner 
so that the chains may be described by a sequences of 
such spacings M, Mi, M2, . . . , Mp. We restrict the form 
of the randomness to cases in which the same sequence 
can be used to describe both chains. The advantage of 
this restriction is that the interchain interaction part of 
the Hamiltonian is still diagonal in the basis of chain 
symmetrized / antisymmetrized wavefunctions . 

When considering the chain- 

symmetrized/antisymmetrized subspaces individually, 
the Hamiltonian given by Eqs. I2I3I is identical to the 
two-component random lattice problem, which has been 
extensively studied [2^ 0, |2^ [l^ . For completeness we 
recapitulate the discussion of the general formalism of 
the coherent potential approximation (CPA) developed 
therein while applying that formalism to the system at 
hand. We begin by expanding the propagator of the 
full, disordered system G in terms of the propagator of a 
uniform system Go with a fixed, and as yet unspecified 
on-site potential. The propagator for the uniform system 
is given by 

Go = (29) 



E-H. 



where 



N 



N 



H, = -tY^{\i+l){i\ + \t){(i+l\)+vY,\t){£\ (30) 



1=1 
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and V is the spatially independent on-site energy. The 
full Hamiltonian of the disordered system is given by the 
sum of Hq from Eq. 1301 and scattering potentials at each 
site given by 



N 



(31) 



In the above equation, the random variable introduces 
the quenched disorder by introducing uncorrelated tun- 
nelling sites with density p along the chain. This random 
variable is takes the value A with probability p and is 
zero otherwise where, as before, A is ±t' for the chain- 
antisymmetrized / symmetrized Hamiltonian. 

One can show in the usual way that the propagator 
for the full, random system G can be expressed pertur- 
batively in terms of the tunnelling Hamiltonian Hj and 
propagator of the uniform system so that 



G 



1 



E-H 



Go + GqHjGq 



(32) 



Our goal, as before is to compute the density of states of 
the system, g(E), which is given in terms of the propa- 
gator as 



g{E) = --ImY^me) 



(33) 
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E/2t 



-0.006 



-0.01 




0.2 0.3 

density 



FIG. 5: Plot of the change in electron energy as a function 
of linker density for cases of quenched, disordered positions 
(short dashes) and ordered positions (solid line). By assum- 
ing the independence of the binding sites so that each linker 
contributes Ai^ — {At^ + t''^y^^, one arrives at a reasonable 
linear approximation to the binding energy as a function of 
cross- link density (long dashes). Inset: Plot of the change in 
energy for an ordered lattice vs. 1/density. Both plots use 
X/2t = 0.25. 



where the sum is over a complete set of states. To approx- 
imate this sum accurately, we follow the work of Soven 
|26| and first reorganize the perturbation series given in 
Eq. 021 by introducing the T operator representing the 
Born-series for scattering off a particular site, £, 



Ti = \e){ee-v){e\ + {si- 
1 

" I- {e\G„\£){ei-v) 
so that Eq. implies 



-vf\£){£\Go\i){i\ 



•(34) 



{e'\G\£) = {e'\Go\i)+Y,i^'\Go\l)Ti{£\Go\l)+ (35) 



We now choose the heretofore arbitrary on-site poten- 
tial. One might imagine that it would be reasonable to 
choose that potential to simply be the mean of the on- 
site energies of the disordered system by taking v = pX, 
however, we are not attempting to average the Hamilto- 
nian over the disorder but rather we will be averaging the 
propagator. Thus, using the CPA approach, we choose 
that on-site potential to enhance the convergence of the 
perturbation series for G given in Eg. 13 51 bv taking v such 
that the average of the T operator over the quenched dis- 
order vanishes. Representing the disorder averages by [•], 
we choose v such that 



+ Yl}^'\Go\em£\Go\£)T-^{!\Go\i) + ■■■. 

1=^1 

The successive terms in the above equation represent the 
propagator of the electron in the uniform system followed 
by a correction due to the interaction of that electron 
with the scattering potential at site ^, followed in turn 
by the interaction of that electron with sites £ and £\ the 
higher order terms (not shown above) have an analogous 
interpretation. It is important to note that, due to the 
reorganization the perturbation series in terms of the T 
operators, the effect of the electron interacting with the 
same site successively has been already taken into ac- 
count. Thus, in the sums representing the second and 
higher order terms it is necessary to restrict the summa- 
tion to avoid revisiting the same site twice in a row, e.g. 
£ ^ £ in the second order term above. 



l + v(£\Go\£) "l-(A-t-)(£|Go|£) 



(36) 

Now, the reorganized perturbation series in Eq. 1351 is 
greatly simplified. Since [(£|T^|^)] = and T operators 
at different sites are uncorrelated, [{£\Ti\£){£'\Te,\e')] = 
[{£\Ti\£)] ■ [{£'\Te'\£')] = 0, we see that first correction to 
Go is fourth order in the T operators. In the limit of 
weak scattering, [G] ~ Go (with the advantageous choice 
for V discussed above) makes an excellent approximation 
since the first correction is of order (pA/(2t))'*. 

Finding the density of states in the disordered sys- 
tem requires that we evaluate the trace of the disorder- 
averaged propagator (Eq. I33|) as a function of v; we do 
this in the position representation by writing propagator 
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in the uniform system as 



■K / a 



^ika{l — n) 



{£\Go\m) ^Go(i,m;E) ^ — dk- „ , 

^ ' ' ' >JV > 1 ; 2n J^/^ E-v + 2tcos{ka) 

(37) 

and extracting the diagonal element 



-1/2 



(38) 



and the sign above is determined by whether the magni- 
tude of {E — v)/2t is greater or smaller than unity. In 
the uniform system the sum on position eigenstates is 
trivial as seen by the ^-independence of the RHS of the 
above equation. However, determining Go{£,£;E) and 
V simultaneously requires that we solve two polynomial 
equations given by Eqs. EElandEHl 

The binding energy is determined by integrating Eg. 1331 
over the filled states. For convenience, we chose the Fermi 
level to again lie at Ep = 0. The results are shown in 
Figs. El As with the ordered lattice of binding sites, 
we find that the binding energy of the disordered braided 
chains is very nearly a linear function of the density of 
sites. 

We can conclude that the primary effect of introduc- 
ing binding points between the polymers is the creation 
of "impurity" bands centered at the energy given by Eq. 
El The states in the impurity bands are created at the 
expense of states near the band edges, E = ±2t. How- 
ever, the density of states far from the band edges is 
essentially unchanged. This can be easily proved for a 
random arrangement of the impurities. We first note 
that the density of states is given by 



9iE) = - — 



1 9$ 
ndE 
1 dk 
^dE 



N + Np 



de 

dk 



(39) 
(40) 



where $ is the phase of the wavefunction at the last site 
on the chain, p is the density of impurities, and 9 is the 
average phase-shift as the wavefunction passes through 
an impurity. For an impurity A at the origin {£ = 0), the 
wavefunction takes the form 

|fe) ^^cos{kla + (j))\l)+A^cos{kla + (l) + e)\l), (41) 



where A is an amplitude to be determined. The phase- 
shift is given by 



tan ^ I tan(0) 



A 



t sin(fca) 



- 0. (42) 



Averaging over the incident phase, (j), we find the average 
phase-shift 



tXcos{ka) 
4t2 sin2(fca) + A2 



(43) 



This shift is odd in the impurity strength. A, therefore 
the change in the density of states from the symmetric 
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FIG. 6; Plot of the change in electron energy as a func- 
tion of the interaction strength for random disorder (short 
dashes), perfectly ordered potential (solid line), and under 
the independent-linker assumption where each cross-link con- 
tributes At^ - (At^ +t''^y^^ (long dashes). Here the cross-link 
density is 0.1. 



and anti-symmetric bands will exactly cancel. Therefore, 
the reorganization of electronic energy must occur where 
this argument breaks down, namely when the separation 
between binding sites is comparable to the wavelength, 
i.e. when sin(/ca) ~ 0. At this region near the band edges 
states are removed to form the imaginary wavenumber 
impurity states. 

The role of the impurity band as the source of the 
binding energy is further supported by the "single state" 
or independent binding site approximation (long dashes) 
line in figures \S[ and j^l This line is the binding energy 
that would be obtained if the sole effect of adding a bind- 
ing site were to move a single electron from the bottom of 
the band {E = —2t) to the energy of the isolated bound 
state {Eb = -(4t^ -t- t'^Y^^). This Hne is an excellent ap- 
proximation for the binding energy of both the ordered 
lattice of binding sites (solid line) and the random ar- 
ray of binding sites (short dashes) which supports the 
suggestion that the creation of the impurity band is the 
dominant effect in the reorganization of electronic energy. 



D. Parallel Configurations 

If the chains adopt a parallel configuration, then there 
will be an overlap between each site and the correspond- 
ing site on the opposite chain. This interaction is eas- 
ily diagonalized with the transformation given by Eq. 01 
This yields the dispersion relation 



E±{k) ^ -2t cos k±t'. 



(44) 



We can solve for the binding energy of the polymers 
in this zipped configuration if the polymers are lightly 
doped (a general solution for the zipped configuration 
eigenstates may be found in Appendix A) . If each poly- 
mer has (5/2 electrons and S <^ N then, prior to zipping, 
all electrons will have energy very near —2t. If the poly- 
mers bind over a length of m monomers, then states are 
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created with energy E < —2t with a density of states 
m(4i2 - {E + f )2)~i/^/7r. The electrons in these states 
have lowered their energy by 



rn 

TT 
TO 



AE = - — y/W~pItTW 



{2t - t') 



2t 



(45) 



This corresponds to a binding energy of ~ — 2i''^/^/37r\/t 
per site. 

If the chains zip past a critical length TOc = 
Stt/ cos~^{t' /2t — 1) then all the electrons have energy 
< —2t. These electrons now have imaginary wavenum- 
bers in the unzipped region of the chain, therefore the 
unzipped regions will no longer be metallic. The binding 
energy is now 



. _ 2tm ,6t:^ ^ , 

AE = sin(— - 6t' 

TT m 



2tS. 



(46) 



Since the derivative of this binding energy with respect to 
the length of the binding region (to) is negative, we note 
that the electrons exert an effective pressure to increase 
the size of the zipped region When they form, we 
expect the length of the zipped regions to be limited only 
by the number of conduction electrons. 

The binding energy per binding site in a zipped region 
is ^ i'^^/'^t^^/^ whereas in the braided configurations the 
binding sites reduce the energy by ~ t'^t~^. The zipped 
configuration exhibits stronger binding per site by \ftjt' ; 
since generically t' < t, the electronic degrees of freedom 
favor zipped configurations of the chains over braided 
ones. To determine whether a pair of chains will, how- 
ever, form such zipped regions or the random braids, we 
must consider the energetics associated with the confor- 
mational degrees of freedom of these charged polymers 
in water. 



III. CONFORMATIONAL DEGREES OF 
FREEDOM OF THE CHAINS AND 
EQUILIBRIUM STRUCTURES 

The equilibrium state of aggregating polymers depends 
on the combination of the binding energy due to the elec- 
tronic degrees of freedom and the change in the transla- 
tional/conformation free energy of the polymer. From 
the interplay of these two contributions to the total free 
energy, we explore the transition from free chains in so- 
lution to chain aggregates and discuss their equilibrium 
structure. We compare the free energies of (i) isolated 
chains in solution, (ii) braided pairs of chains consist- 
ing of unbound loops between isolated crossing points, 
and (iii) parallel or "zipped" configurations in which 
the chains have numerous consecutive tunnelling sites 
along their length. The aggregation of DNA in the pres- 
ence of divalent ions has been approached in a similar 
manner p^. however, in the present system there is no 



quantity analogous to the fixed concentration of linking 
divalent ions since the interchain binding sites are created 
spontaneously at chain intersections. 

We begin by considering the coexistence of isolated 
chains and paired ones in solution. The thermodynamic 
potential of a solution with a concentration of chains n 
can be written up to a trivial constant as 

E{Tt2) 

— — — = (n-2n2)(ln (n - 2n2)-l)+n2(lnn2-l)+ri2 
Kb J 

in terms of the concentration of bound-chain pairs n2 [2 
In the above expression 5*0 is the chain configurational 
entropy of one polymer: So cx ksL/lp. Its appear- 
ance in Eq. 07| comes from the fact that in order to 
pair two chains, one of them must adopt the same con- 
figuration as the other and thereby lose its conforma- 
tional entropy. The first two terms of Eq. 1471 represent 
only the translational entropy of the isolated chains and 
bound pairs; they do not incorporate this effect. Min- 
imizing Eq. 1471 we find the concentration of dimers is 
n^exp((— rS'o — £b)/kBT). It remains now to calculate 
the interchain binding energy £b, which is a combination 
of electrostatic chain repulsion and chain binding due to 
the proposed tunnelling mechanism. We do not include 
van der Waals interactions between chains. This longer- 
ranged attractive interaction is generically stronger for 
the conducting polymers since the extended electronic 
states enhance their low-frequency polarizability. At 
higher frequencies one expects that the polarizability at 
higher frequencies is typical of small organic molecules 
since the dielectric spectrum at higher frequencies is pri- 
marily due to localized molecular states that are similar 
to all hydrocarbons. Such van der Waals interactions 
lead to a longer range attraction that, in high salt con- 
centrations at least, are a subdominant correction to the 
interchain binding potential when the interchain distance 
is on the order of Angstroms. In addition the enhance- 
ment of the Van der Waals interaction is significant for 
only the metallic rather than the semiconducting state of 
the polymers. Thus we expect that the interchain tun- 
nelling mechanism discussed in this article to dominate 
the binding energy of the observed bundles. We examine 
this point further in appendix B. 

We first consider the binding energy of the more loosely 
bound, braided configuration of the molecules, £^°-'^<^ , 
The binding energy is a combination of the electrostatic 
interaction of the chains, the bending of the molecules on 
scales short compared to their persistence length, and the 
electronic binding energy of the crossing point. At high 
salt the first of these terms is significant only near the 
crossing points. There the charged chains can be treated 
as crossed straight rods. If two charged rods cross at an 
angle Q with a distance of closest approach d the electro- 
static energy is 

E,, . H^Sg^l. (48) 




Kat. sm f 



where Ib 



^ /eksT is the Bjerrum length, Oc is the 
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average separation between charges along the backbone, 
and is the Debye screening length. 



The bending energy of a polymer is 



Ebe 



ksTL 



dl 



my 



(49) 



where Ip is the persistence length, and R{1) is the local 
radius of curvature. If the separation between interaction 
points is comparable to the persistence length, then in 
the limit Ip the crossing angle will be very close 

to tt/2 and the arc the polymers follow will be very nearly 
circular. Then the energy per length of the polymers, L, 
is 



Ebraid{l)/L — ( 



-Kd fl2 hoTl TT^ 



2'KkBTlBe 



where / is the arc length of the polymer between intersec- 
tions. Minimizing with respect to I, the energy becomes 



^braid 



"iTTkBTlBer 



2t I 



-L. 



(51) 



Here we have neglected the configurational entropy of the 
loops of chain making up the braid. As long as I < Ip we 
expect this entropic contribution to be negligible. 

We now turn to the zipped configuration of the chains. 
To consider the electrostatic interaction between the two 
chains, we note that the potential, <i>, around a thin 
charged rod in a salt solution is given by 



2lBkBT 
$(r) = Ko{Kr), 



(52) 



so the energy of two polymers separated by a distance d 
over a length L is 



2lBkBT / ,N r 2t . 5nat 
L 5 Ko(Kd)- 



-L sm , 



-t'S+2tS 
(53) 

where Kq is the modified Bessel function, the distance 
between tight-binding sites is at, and we have assumed 
the binding energy is given by Eq. 1461 

We observe that both Eqs. 1511 and 1^ depend on the 
length over which the binding occurs. So, before we use 
either of these expressions in Eq. E|we must choose the 
suitable length, L*. For braided chains L* is the length 
of the entire chain and for zipped chains L* is determined 
from E3 yielding L* ~ 2tdal/lBkBTKQ{Kd). Note that 
^0 in Eq. ^\is also proportional to L* . 

Based on these considerations we propose a schematic 
phase diagram for charged, conducting polymers in solu- 
tion. This diagram (shown in figure [TJ, which represents 
the range of phase behavior of these polymers, is spanned 
by the persistence length and the (suitably scaled) in- 
verse, linear charge density of the chains. This inverse 
charge density is measured in terms of the ratio of the 
distance Uc between charged groups along the polymer 
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(a.u.) 
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FIG. 7: Schematic phase diagram for two charged, conduct- 
ing polymers. This diagram is spanned on the vertical axis 
by the persistence length scaled by the total contour length 
of the chain while the horizontal axis is the dimensionless ra- 
tio of the the distance between charges along the chain Oc 
to the distance between potential tunnelling sites at. In the 
upper left region of the figure, chains of high charge density 
and long persistence length are unbound in solution (free); 
upon reducing the linear charge density along the chain and 
moving to the right in the diagram the chains generically be- 
come bound in a parallel or zipped configuration. For flexible 
enough chains of intermediate charge density, however, more 
loose, braided configurations are predicted. 



backbone divided by the distance between potential tun- 
nelling sites, at- In the upper right portion of the figure 
the chains are both stiff and have a high charge den- 
sity relative to their density of potential binding sites. 
Thus, they are unbound in equilibrium. Generally, mov- 
ing down and to the right in the diagram induces inter- 
molecular binding due the mechanism presented in this 
work. The equilibrium bound configurations take the 
form of parallel, tightly bound pairs, which we denote as 
zipped chains or as more loosely bound braids. 

The onset of braiding is determined by setting that 
contour length between binding sites equal to the to- 
tal contour length of the chain, L. This gives the free- 
chains/braid boundary as 



L 



2t 



2T:k^Tlpe 



— Kd 



(54) 



The term in parenthesis on the right-hand side of the 
above equation is collectively the effective binding of one 
interchain crossing. The first part of that term is the 
binding energy due to the reorganization of the electronic 
states of the polymers while the second term represents 
the reduction in binding energy associated with the elec- 
trostatic interaction of the two chains. For the cases of 
current interest in which a single interchain bond is fa- 
vorable the right-hand side of Eq. |^is positive. 

To explore the boundary between zipped and braided 
states of the polymers we compare the energy of the 
braided configuration given by Eq. 1511 to the energy of 
the zipped configuration. This latter energy is a combi- 
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nation of the binding energy of a zipped region given by 
Eq. 1^ For the case of hghtly doped chains so that con- 
tour length of the zipped regions L* is smaUer than that 
of the whole polymer, one finds that the boundary be- 
tween braided and zipped chain configurations to occur 
where 



L 



,T7r2 ( 2t'3/^ 



2lp fee TKq {Kd)ac 



(55) 



This phase boundary continues to smaller ad at until 
it intersects the free/braid boundary discussed above. 
There it terminates. Finally, there is a vertical line sepa- 
rating free and zipped chains that intercepts the intersec- 
tion point of the free/braid and braid/zipped boundaries. 
This last separation is independent of chain persistence 
length since the energy of neither the zipped configu- 
rations nor the free chains depends on the persistence 
length. 



IV. DISCUSSION 

We find that conducting polymers in the metallic state, 
i.e. where the Fermi level lies within a band, will gener- 
ically aggregate in solution due to the reorganization 
of the electronic degrees of freedom of these molecules 
upon the close (Angstrom-scale) approach to each other. 
This reorganization of the extended electronic states of 
the molecule results in the formation of localized bound 
states at these crossing-points and the consequent reduc- 
tion of the electronic energy more that compensates for 
the reduction of chain configurational entropy due to the 
creation of the cross-links between the two polymers. It 
is interesting to note that these localized bound states are 
the direct analogs of the standard binding/ anti-binding 
orbitals that are created out of atomic electronic states 
upon the close approach of their respective nuclei. In 
the case of conducting polymers, however, these local- 
ized bound states are created out of the extended states 
of the metallic molecules. 

The localized nature of these states has the conse- 
quence that the binding energy is very nearly linear in 
the density of binding sites when these sites are separated 
by at least a few monomers as they are in all braided 
configurations. However, when the polymers lie in par- 
allel, zipped configurations, the resulting binding energy 
per site is stronger than the isolated binding sites by 
0{{t/t'Y/'^). The equilibrium structures formed by these 
polymers will depend not only on the interchain binding 
mechanism, but also the linear charge density of the poly- 
mers, solvent screening, and the persistence length of the 
polymer. The combination of these factors will determine 
whether the polymers form loosely braided structures, or 
tight bundles. 

Based on these considerations we have proposed a 
schematic phase diagram showing the expected equilib- 
rium structures that may be observed in such systems. 



For large enough linear charge density, one finds that the 
lowest energy state of the aggregate is a braid consist- 
ing of loops of unbound polymer between cross-linking 
binding sites. For smaller linear charge density one finds 
that the energetically favorable configuration is two par- 
allel chains where the interchain distance remains on the 
order of a few Angstroms all along their arc length. 

Clearly, the direct experimental observation of such 
structures is difficult. We note that SANS experiments 
llTj are consistent with the formation of bundles in 
aqueous solutions, but the spatial resolution of these ex- 
periments is not fine enough to distinguish between the 
braid and parallel phases discussed above. Moreover, the 
current calculations are not directly applicable to these 
experiments since previous experimental studies focussed 
on semiconducting polymers. Nevertheless, our proposed 
binding mechanism mediated by the reorganization of the 
electronic degrees of freedom is also relevant to the semi- 
conducting casej^. Since it will likely remain difficult 
to quantitatively test the detailed structure of the ag- 
gregate one must directly probe the electronic states by 
studying the spectroscopic signature of the molecules as 
they aggregate. 
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APPENDIX A: ZIPPED CONFIGURATION 
STATES 



Here we solve for the eigenstates for two chains that 
are zipped together over a portion of their length. The 
Hamiltonian is 



N-l 



N-1 



H, = -tY^{\l+l){t\ + \l){i+l\)+Y^ vM){t\ (Al) 



1=1 



1=1 



where = in the unzipped regions (1 < / < n and 
n + m + 1 < I < N) and tj; = A in the zipped region 
(n -f 1 < ^ < n + to). Here A — ^t' for the symmetric 
(anti-symmetric) states, to is the number of sites in the 
zipped region, and n{n') is the number of sites in the 
unzipped region to the left(right) of the zipped region. 
The chains have a total length of = n -I- to -I- n'. The 
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eigenfunctions are of the form 



|fc) = ^sin(Zfc)|0+A sm{lk' + (j))\l) 

1=1 l=n+l 
N 

+B J2 sm{lk + e)\l), (A2) 

l=n-\-Tn-\-l 

where A and B are undetermined amphtudes, and 
— 2icosfc = —2tcosk' + A. The boundary conditions are 

sin((n+l)A;) = Asm{{n + l)k' + (/)) 

sin(nfc) = Asiii{nk' + (f)) 

Asm{{n + m + l)k' + (/)) = B sm{{n + m + l)k + 6) 

A siii{{n + m)k' + (f>) = B sm{{n + m)k + 6) 

{N+l)k + e ^ jn, (A3) 

where j is an integer. After some algebra, we find two 
equations for k and 4> 



sin(fcn') 



sin(fc'(n + m + 1) + 



sin(fc(n' + 1)) sin(A:'(n + m) + 0) 

cot(fcn) sin(fc) — A 



cot(A:'n + cj)) 



sin(A:') 



(A4) 
(A5) 



From these equations we find the quantization condition 
for k 

CC' + (? + C) sin(fc') cot(fc'(m + 1)) - sm^{k') = 0, (A6) 
where 



^ = sin(A;) cot(nfc) — A 
^' EE sin(A;) cot(n'fc) - A. 



(A7) 



From these conditions (Eqs. IA6I IA7(I and the dispersion 
relation 



Ek = — 2icos(fc), 
one obtains the energy eigenvalues. 



(A8) 



APPENDIX B: INTERMOLECULAR VAN DER 
WAALS ATTRACTION 

The van der Waals interaction energy between two like 
molecules in a screened medium is given by 



2 -2«;r 



(Bl) 



where ao is the polarizability of the molecules, e is the 
dielectric constant of the medium, is the Debye 

screening length, and / is the ionization energy of the 
molecules '31]. At wide separations the polymers can be 
approximated as conducting spheres with a radius 
thus ao — R^g- In aqueous solutions e ~ 80 and for an 
/ of order eV then Ei,dw ^ leV since r ^ Rg. If, on 
the other hand, r < Rg it is more accurate to describe 
the polymers as narrow ellipsoids. In this situation we 
approximate the polymers as ellipsoids of length Ip and 
thickness a that is a monomeric dimension. The polariz- 
ability along the wire is found to be of order [s^- 



" In^ 



(B2) 



and perpendicular to the wire it is 



27 



(B3) 



So, the interaction energy for parallel and perpendicular 
wires is 



En 



E, 



76 -2Kr 



'- P 



(B4) 
(B5) 



For Ip ~ lOnm, e ~ 80, a ~ r ~ l/2nm, / ~ leV, 
and ~ Inm, we find that E^^ ~ lO^'fcsT and E± ~ 
2kBT. At high salt concentrations, when the electro- 
static screening length k"'^ ~ 2 A, these dipolar interac- 
tions are suppressed by a factor of ~ 10** suggesting that 
the binding mechanism discussed in the article dominates 
at short interchain separation. At lower salt concentra- 
tions van der Waals interactions play a larger role in the 
formation of the aggregates, but the structure of those 
aggregates will still depend on the interchain electron 
tunneling mechanism. 
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